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Abstract 

We consider discrete spin models on arbitrary planar graphs and 
lattices with frustrated interactions. We first analyze the Ising model 
O | with frustrated plaquettes. We use an algebraic approach to derive 

the result that an Ising model with some of its plaquettes frustrated 
CN ■ has a dual which is an Ising model with an external field iir/2 applied 

| to the dual sites centered at frustrated plaquettes. In the case that 

■ all plaquettes are frustrated, this leads to the known result that the 

dual model has a uniform field m/2 whose partition function can be 
evaluated in the thermodynamic limit for regular lattices. The analysis 
{Sj ' is extended to a Potts spin glass with analogous results obtained. 
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1 The frustrated Ising model 



A central problem in the study of lattice-statistical problems is the consid- 
eration of frustrated spin systems (see, for example, [1] - [4]). A particularly 
useful tool in the study of spin systems is the consideration of duality rela- 
tions (see, for example, [5, 6]). Here we apply the duality consideration to 
frustrated discrete spin systems. 

We consider first the Ising model on an arbitrary planar graph G. A 
planar graph is a collection of vertices and (noncrossing) edges. Place Ising 
spins at vertices of G which interact with competing interactions along the 
edges. Denote the interaction between sites i and j by — = —SijJ where 
Sij = ±1 and J > 0. Then the Hamiltonian is 

H(a;S) = - SijJ (!) 

<i,3> 

where <7j = ±1 is the spin at the site % and the summation is taken over all 
interacting pairs. 

The Hamiltonian (1) plays an important role in condensed matter physics 
and related topics. Regarding as a quenched random variable governed 
by a probability distribution, the Hamiltonian (1) leads to the Edwards- 
Anderson model of a spin glasses [7]. By taking a different SV,-, the Hamilto- 
nian becomes the Hopfield model of neural networks [8]. Here, we consider 
the Hamiltonian (1) with fixed plaquette frustrations. 

Let G have N sites and E edges. Then it has 

N* = E + 2-N (Euler relation) (2) 

faces, including one infinite face containing the infinite region and N* — 1 
internal faces which we refer to as plaquettes. The parity of a face is the 
product of the edge factors around the face which can be either +1 or 
— 1. A face is frustrated if its parity is —1. An Ising model is frustrated if 
any of its plaquettes is frustrated, and is fully (totally) frustrated if every 
plaquette is frustrated. The fully frustrated model is also known as the odd 
model of the spin glass [1]. As the parity of the infinite face is the product 
of the parities of all plaquettes, the parity of the infinite face in a totally 
frustrated Ising model is —1 for N* = even and +1 for N* = odd. An 
example of a full frustration is the triangular model with = —1 for all 
nearest neighbor sites 
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The values of parity associated with all plaquettes defines a "parity con- 
figuration" which we denote by T. The set of interactions {SV,} corresponding 
to a given T is not unique. For the triangular model, for example, any {SV,} 
which has either one or three = — 1 edges around every plaquette is totally 
frustrated. For a given {SV,} and T, the partition function is the summation 

z{{Sij})= E ••• E IK-"''" (3) 

(71=±1 (7JV=±1 E 

where the product is taken over the E edges of G. 
1.1 Gauge transformation 

A gauge transformation is site-dependent redefinition of the up/down spin 
directions. Mathematically a gauge transformation transforms the spin vari- 
ables according to [2] 

<?i -> (J- = Wi (Ti, i = 1, .., N. (4) 

In the above if Wi = +1 the original definition of up/down spin directions is 
maintained, and if — — 1 the definitions of up/down are exchanged. Under 
the gauge transformation the in Eq. (1) transform as follow 

Sij -> S^- = Wi S i:j wj V (5) 

Since w 4 2 = 1, we have 

H{a;S) =H{a'-S'). (6) 

Clearly, the gauge transformation (5) leaves the parity configuration Y un- 
changed, i.e., 

II ^ IK Vface. (7) 

face face 

For each parity configuration V there are 2 N ~ l different {Sij} patterns 
consistent with it. To see that we note in Eq. (5) each of the 2 N choices of 
{wi} leads to a new {S^} except the negation of all Wi which leaves {SV,} 
unchanged. Conversely, any two sets of interactions {5V,} and {5^} for the 
same T are related by a gauge transformation which can be constructed as 
follows. Starting from any spin, say, spin 1, assign the value w± = +1. One 
next builds up the graph by adding one site (and one edge) at a time. To the 
site 2 connected to 1 by the edge {12}, one assigns the factor w 2 = W\ Si 2 S' 12 
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which yields Wi S i2 w 2 = S[ 2 consistent to (5). Proceeding in this way around 
a plaquette until an edge, say {nl}, completes a plaquette. At this point one 
has 

w n s nlWl = ( n n s«W = ^ (8) 

\ plaquette / \ plaquette / 

which is again consistent to (5). Continuing in this way one constructs all u>j 
which transform {S^} into {S^}. Note that if we had started with w\ = — 1, 
we would have resulted in the negation of all Wi. Thus, the bijection between 
the 2 N ~ l sets {5V,-} and 2 JV_1 gauge transformations is one-one. 

In addition to Eq. (7), the gauge transformation also leaves the partition 
function invariant [9, 10], i.e. 

z{{s i3 }) = e n = e n eS ' ijJ<a ' j = z ({^}) • ( 9 ) 

{a} E {a'} E 

As a result the partition function only depends on Y and we can rewrite (3) 

as 

z(r) = 2 -("-U (10) 

{Sij} 

where the summation is over all 2 JV_1 distinct {Sij} consistent with the par- 
ity configuration T for the same partition function. This expression of the 
partition function is used to derive the duality relation in ensuing sections. 

1.2 The fully frustrated Ising model 

For our purposes it is instructive to consider first the case of full frustration. 
Duality properties of fully frustrated model have previously been considered 
by a number of authors [2, 4] for regular lattices. We present here an alternate 
formulation applicable to arbitrary graphs and arbitrary frustration. 

The graph D dual to G has iV* sites each residing in a face of G, and 
E edges each intersecting an edge of G. We restrict to iV* = even so that 
all faces of G including the infinite face are frustrated. This restriction has 
no effect on the taking of the thermodynamic limit in the case of regular 
lattices. Since the signs around each face are subject to the constraint 
n = —1, we introduce in the summand of (10) a face factor (1 — 11 Sij)/ 2 
and sum over S^ = ±1 independently. Similarly, writing Oij = di&j, we can 
replace summations over cTj = ±1 in (10) by = ±1 by introducing a factor 
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(1 + H<Tij)/2 to each face. Thus, (10) becomes 
- 2 . 2 - ( iv-i) 



' FF 



e E[n^ j ]n[a+n^)(i-n^)] (n) 

Wij} {Sij} E face 



where the subscript FF denotes full frustration, and the extra factor 2 in (11) 
is due to the 2 — > 1 mapping from ai to a^. 

For a face having n sides we rewrite the face factors as 

i + LK = EI1%^) (12) 

i-n^ = eii^^) as) 

where each product has n factors 

F(a; fx) = <5 M+ + a 5^ 

G(S; u) = 5 V+ + S uo n S v - , (14) 

5 is the Kronecker delta function and cu n = (— l) -1 /™ = e _l7r//n . 

We now regard ji and v as indices of two Ising spins residing at each dual 
lattice site. After carrying out summations over and S^, the partition 
function (11) becomes 

^^^TEn^ft^v') as) 

where we have made use of the Euler relation (2) and B is a Boltzmann 
factor given by 

B(jm,u;im',i^= E E e s ° J F(a;Li)F(a;Li')G(S;v)G(S;v'). (16) 

a=±l S=±l 

Here G(S; v') is given by (14) with uj n — > u n > = e~ m l n ' and the two faces 
containing spins {//, z/} and {//, z/} have respectively n and n' sides. 
Substituting (14) into (16) and making use of the identities 

<Vh<V+ + <V-<V- = (1 + /V)/ 2 

<WV- + W+ = (i-/V)A (17) 

one obtains 

£(//, v\ li', v 1 ) = 2A{1 + lili') cosh J + 25(1 - lili') sinh J (18) 
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where 



B = u n 5 v -5 u i + + uj n i5 u+ 5 l/ i 



(19) 



Number the four states {//, is} = {+, +}, { — , — }, { — , +}, {+, — } respectively 
by 1, 2, 3, 4. The Boltzmann factor (18) can be conveniently written as a 4 x 4 
matrix 

/B n B 12 \ 
B21 B22 
B n B l2 
V B 21 B 22 J 

where 



B(n,v;n',i/) 



(20) 



Bn = 4 cosh J, i?i 2 = 4 av sinh J, 

B 2 i = 4 w n sinh J, £> 2 2 =4w n cosh J . (21) 

Thus the partition function of the {/i, u} spin model is twice that of an 
Ising model on the dual lattice. The exchange coupling constant K and the 
magnetic field h in the dual model are determined by 

B U = De K+{h/n)+{k ' /n '\ B 12 = JJ e - K Hh/n)-(h'/n') 

B 21 = De -K-(h/n) + (h'/n')^ ^ = De K-(h/n)-(h>/n>) _ ^ 

Here n and n! are the number of edges incident at the two dual sites respec- 
tively. 

The solution of the above equations gives 

e~ 2K = tanhJ> 0, D = 4(u n ccV^Vsinh J cosh J 

e 2(h/n) = l/Un = ^/n ^(fc'/n') = = ^/n> ? (33) 

or equivalently 

I / \ 7 7T 

K = — ln(tanhj) and h = h' = —. (24) 
Thus, we have established the equivalence 

Z FF (J) = 2 7V " 1 i N * (sinh J cosh J) E ' 2 if) (25) 
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where Zj sin g(i7r/2, KJ is the partition function of a ferromagnetic Ising model 
on D with interactions K > and an external field in/2. In writing down 
(25) we have made use of the identity 2 • 2~ ( - E+N *^ A E = 2 N ~~ 1 and the fact 
that (uj n uo nl ) E ' 2 = (-i) N * = i N * for N* = even. 
Remarks: 

1. The duality relation (25) has previously been obtained by Fradkin et 
al. [2], and for the square lattice by Suzuki [4] and Siito [11], and by Au-Yang 
and Perk [12] in another context. 

2. The duality relation (25) is different from the Kadanoff-Ceva-Merlini 
scheme [13, 14] of replacing K by K + in/ 2 (corresponding to J < in 
(23)) in the ferromagnetic model. Suzuki [4] has made the explicit use of the 
Kadanoff-Ceva-Merlini scheme in deriving (25) for the square lattice. For 
fully frustrated systems the Suzuki method can be extended to any graph 
whose dual admits dimer coverings. 

3. The duality relation (25) holds for a fixed {SV,} without probability 
considerations, and therefore differs intrinsically from that of a spin glass 
obtained recently by Nishimori and Nemoto [15] using a replica formulation. 

4. The duality relation (25) which holds for any lattice appears to support 
the suggestion [3] that all fully frustrated Ising models belong to the same 
universality class. 

1.3 The thermodynamic limit 

The partition function (25) for an Ising model in a uniform field in/2 can be 
exactly evaluated for regular lattices. Defining the per-site "free energy" 

Lee and Yang [16] have obtained a closed form expression of f(K) for the 
square lattice. Their result, which was later derived rigorously by McCoy 
and Wu [17] and others [14, 18], is 

/ = i- + C + -^ r d6 r d<f)\n[z + z' 1 +2-4cos0cos0], (27) 

2 167T 2 J-7T J-W 

where C = [ln(sinh2K)]/2, z = e~ AK . The free energy (27), which is the 
same as that obtained by Villian [1], is analytic at all nonzero temperatures. 

The solution for the triangular Ising model in a field in/ 2 has also been 
deduced previously [18, 19]. However, it can also be obtained most simply 
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by observing that the honeycomb lattice, which is the dual of the triangular 
lattice, has a coordination number 3. It follows that we can recast the field 
Boltzmann weights as e ma ^ 2 = io~j = io~j and redistribute the cr| factor at 
site j to its three incident edges. Then, as pointed out by Suzuki [4], each 
edge can be associated with a factor iaiaje K = e ( K + m / 2 )< J i< J i anc [ the desired 
solution can be obtained from that of the zero-field honeycomb lattice with 
the simple replacement K — » K + in/2. This gives 

f = i^ + C+ — ^ r old r d(p In [(1 + e 4K ) 2 + 4 cos 0(cos 9 + cos 0)1 , (28) 

Z 107T J— 7T J— 7T 

where C = [ln(2 sinh2i^)]/2 and K is the Ising interaction on the honeycomb 
lattice. Again, there is no finite temperature phase transition. 



1.4 Arbitrary plaquette parities 

In a similar fashion one can extend the above analysis to Ising models with ar- 
bitrary face parities. All steps of previous subsections can be carried through 
except that for faces which are not frustrated we must replace u) n by 1 at the 
corresponding dual sites. This results in a zero field (instead of a field in/ 2) 
at these sites. Thus, for an Ising model with arbitrary parity configuration 
T, its dual model has fields and in/2, respectively, at sites in faces of parity 
+ 1 and —1. Explicitly, we find 

Z{V) = 2 iV - 1 H)^(sinh Jcosh J) £/2 4m ) g ({^},^) (29) 
where the dual partition function is 

ziSl g {{hj},K) = J2U^ j U^ m ■ (30) 

{fj.i} E face 

Here, Np is the number of frustrated faces and the external field at site j 
is hj = in/ 2 or depending on whether the face associated with the site is 
frustrated or not. 
Remarks. 

1. The duality relation (30) for Ising models with arbitrary frustrated 
plaquettes can be found as contained implicitly in [2]. 

2. By writing e* 7 ™"/ 2 = ia in the dual partition function, we see that the 
partition function of an Ising model with p frustrated faces is dual to a jo- 
spin Ising correlation function in zero field. Particular, the p = 2 correlation 
problem has been studied in details [17] which now leads to a wealth of 
information on the correlation of two frustrated plaquettes. 
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2 Potts spin glass 



Our analysis can be extended to a g-state spin model, the Potts spin glass. 
First, we recall the definition of a chiral Potts model. The chiral Potts model, 
which was first considered in [5] , is a discrete g-state spin model with a cyclic 
Boltzmann factor A(£, £') = A(£ — £') between two spins at sites % and j and 
in states & and ^ = 0, 1, • • • , q — 1. The interactions are g-periodic, namely, 
the Boltzmann factor satisfies 

U(t + q) = U(t). (31) 

The interaction can be symmetric, namely, £/(£) = £/(—£), as in the case of 
the standard Potts model [20], but in our considerations this needs not be 
the case. 

A Potts spin glass is a chiral Potts model with random interactions. To de- 
scribe the randomness one introduces edge variables \j — Xji — 0, 1, • • • , q— 1 
and considers the partition function [15, 21, 22] 

Z Potts ({A,,}) = J2 n Ufa - + Ay) . (32) 

fi=0 E 

Note that if U is symmetric and q = 2, the partition function (32) reduces 
to (3). A plaquette has "flux" r (= 0, 1, 2, • • • , q - 1) if [23] 

Aij = r (mod q). (33) 

plaq 

A set of leads to a flux configuration T, which is specified by the values 
of the flux for all faces. 

2.1 Gauge transformation 

A gauge transformation for the Potts spin glass is the mapping 

-Vj ~~ *■ -Mj = ^ij + 9j, (34) 

where 6i = 0, 1, • • • , q — 1. It is clear that this mapping leaves the flux 
configuration T unchanged, i.e., 

E A« = E ^- (35) 

plaq plaq 
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Since a global change of all £j by the same amount preserves the total 

number of distinct {Ay} consistent with a particular flux configuration T is 
qN-i^ Conversely, any two sets of {Ay} and {A y } giving rise to the same 
flux configuration are related through a gauge transformation. To see this 
we start from a arbitrarily chosen site, say site 1, and set Q\ — 0. Next we 
assign 9 2 = #1 + A i2 — A' 12 to site 2 connected to site 1 by an edge. Continuing 
in this way as in the Ising case, one eventually determines a set of ^ which 
transforms {Ay} into {A y } , and vice versa. The bijection between the g^ -1 
configurations of Ay and gauge transformations for a given T is one-one. 

In addition to leaving the flux configuration unchanged, gauge transfor- 
mation also leaves the partition function invariant, namely, 

-Zpotts ({A i ,})=Z Potts ({A: j }). (36) 

Thus, in analogous to (10), we have 

Z Potts (T) = £ z({Ay}) (37) 

which is used to derive a duality relation. Again, the sum in Eq. (37) run 
through all {Ay} consistent with a given flux configuration T. 

2.2 Duality relation 

In the Potts partition function (32), write £y = £j — £j and to each face 
having a flux r introduce two factors 

1 9-1 

_ y e -i 27 r M (62+63+-+^i)/<? = ! if &2 + • • • + f„i = (mod q) 

= 0, otherwise 

1 9-1 

_ J- e -<^(A 12 +A 2S +-+A f , 1 -r)/, = 1? if Al2 + . . . + Anl = r ( mod g) 

5 i/=o 

= 0, otherwise . (38) 

This permits us to sum over £y and A^ independently. Thus, in analogous 
to (15), we obtain 

^Potts(r) = g- £ -^ £ J2 v\ y') (39) 
M M 
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where 



B(im, v- //, v>) = ]T £ U(£ + A) 



x exp 



i2n ( rv r'v r 
-(n-n)£-{v-v)\ 



n rv 



(40) 



and n, n' are the numbers of sides of the two plaquettes containing {//, u} 
and {//', i/'}, and flux r and r' respectively. 

Carry out the summations in (40) after introducing the Fourier transform 

1 9-1 

U(t + A) = - £ Hv) J 2 ^ +x)/q (41) 

1 r?=0 

where A (77) are the eigenvalues of the matrix U [5]. One obtains 

B(n, is; //, v') = q <W,— ' A(u - v<) e i2 ^ n j^ v >/ q n> . ( 42 ) 

In the above 5^-^^-^ sets [/, — [/,' to v — 1/ (mod q). 

The substitution of (42) into (39) followed by summing over /i now yields 
the result 

Zpotts(r) = g 1 - iY *Z Potts ( D )({^}, A) (43) 

where 

w d) ({m, a) = e [ija^, ^)][n^] • ( 44 ) 

{i/j} _B face 

is the partition function of a chiral Potts model on the dual graph D, which 
generalizes (29) to Potts spin glasses. The dual chiral Potts model has Boltz- 
mann weights A(/ij, jij) — A(//j — /ij) and external fields 

hj = i—^, j = l,2,---,N\ rj = 0,l,...,q-l (45) 

on the spin in plaquette j which has a flux rj. When r,- = for all j, (43) 
reduces to the duality relation for the zero-field chiral Potts model given by 
Eq. (13) in [5]. 
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3 Summary and acknowledgments 



We have obtained duality relations for planar Ising and chiral Potts models on 
arbitrary graphs and with fixed plaquette parity or flux configurations. Our 
main results are the equivalences (25) for the fully frustrated Ising model, 
(29) for the Ising model with arbitrary plaquette parity, and (43) for the 
chiral Potts model with arbitrary flux configurations. In all cases the dual 
models have pure imaginary fields applied to spins in plaquettes that are 
frustrated and/or having a nonzero flux. 

We thank H. Kunz for a conversation which motivated this work. We are 
also indebted to C. Henley, H. Nishimori, J. H. H. Perk, L. Pryadko, and A. 
Suto for calling our attention to relevant references. 

Work has been supported in part by NSF Grants DMR-9971507 and 
DMR-9980440. 
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